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We perform a symmetry analysis of the zero-temperature instabilities of the tetragonal phase of 
Ca3Mn2 07 and Ca3Ti2 07 which is stable at high temperature. We introduce order parameters 
to characterize each of the possible lattice distortions in order to construct a Landau free energy 
which elucidates the proposed group-subgroup relations for structural transitions in these systems. 
We include the coupling between the unstable distortion modes and the macroscopic strain tensor. 
We also analyze the symmetry of the dominantly antiferromagnetic ordering which allows weak 
ferromagnetism. We show that in this phase the weak ferromagnetic moment and the spontaneous 
ferroelectric polarization are coupled, so that rotating one of these ordering by applying an external 
electric or magnetic field one can rotate the other ordering. We discuss the number of different 
domains (including phase domains) which exist in each of the phases and indicate how these may 
be observed. 

PACS numbers: 61.66.Fn,61.50.Ks,75.85.+t 



I. INTRODUCTION 
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FIG. 1: (Color online) The unit cell of the RP (n = 2) system (left) and the (n = 1) system (right). 



Ruddlesden-Popper (RP) systems [l[ are compounds of the form A2+„Bi +n C4+3 n , where n is an integer, and the 
valences of the ions are usually A = +2, B = +4 and C (oxygen) = —2. Systems like SrTiC>3 can be regarded as being 
ri = oo. At high temperatures T rs 800K, their crystal structure is tetragonal, consisting of n-layer units, each layer 
consisting of vertex sharing oxygen octahedra at whose center sit a B ion, as shown in Fig. [TJ (For an alternative 
illustration, see Ref. 0. As the temperature is lowered these systems can undergo structural phase transitions into 
orthorhombic structures. [3j-|7| These structural transitions usually involve reorientation of oxygen octahedra, a subject 
which has a long history, Q of which the most relevant references for the present paper are those (8l-fl0| which deal with 
the general symmetry aspects of these transitions. One reason for the continuing interest in octahedral reorientations 



is because they are important for many interesting electronic properties, such as high- T c superconductivity, [11| colossal 
magnetoresistance. [l2| . metal-insulator transitions. [l3| and magnetic ordering, [lj] 

Here we will focus on the n = 2 RP systems, CasMnoO? (CMO) and Ca3Ti2C>7 (CTO), which present less complex 
scenarios than the n — 1 systems. It seems probable [7U15I] that well above room temperature the crystal structure 
of CMO is that of the tetragonal space group I4/mmm, # 139 (space group numbering is that of Ref. [Lsl ). and 
at room temperature that of space group Cmc2i (# 36). [I?) The theory of isotropy subgroups [l8j strongly forbids 



a direct transition from I4/mmm to Cmc2i. In fact, first principles calculations on these systems by Benedek and 
Fennie (BF) [lj| and on other systems [2(| indicate that the transition from I4/mmm to Cmc2i should proceed via 
an intermediate phase which is probably Cmcm (#63), consistent with Ref. Ila. Up to now no such intermediate 
phase has been observed for CMO or CTO. Unlike the other phases, the Cmc2i phase does not possess a center of 
inversion symmetry and is allowed to have a spontaneous polarization. Recent measurements on ceramic Ca3Mn207 
find a clear pyroelectric signal consistent with the onset of ferroelectric order close to T* = 280K.[2l| Therefore T* is 
identified as the temperature at which Cmc2i appears. Since this ferroelectric transition seems to be a continuous and 
well-developed one and since a direct continuous transition between I4/mmm and Cmc2i is inconsistent with Landau 
theory, [18| the seemingly inescapable conclusion is that the phase for T slightly greater than T* is not I4/mmm, but 
is some other phase which does not allow a spontaneous polarization. Thus the phase at temperature just above 
T = 280K may be the long sought for Cmcm phase. In fact, the Cmcm phase has been observed in the isostructural 
compounds LaCa2Mn 2 07[22] and Bio.44Ca2.56Mn 2 07[23[ at room temperature. In view of the results of Refs. and 
Il5l . it is possible that the Cmcm phase may exist only over a narrow range of temperature. As the temperature is 
further lowered, an antiferromagnetic phase is observed. [UH | In this phase, which appears at T = 115K,[l4j the 
antiferromagnetic order is accompanied by weak ferromagnetism. ^4ll24| 

Theoretically, there have been efforts to understand systems like these from first principles calculations. For 
instance, Ref. [25| found the nearest neighbor exchange J nn within a bilayer to be J nn /fcB = — 39K, giving a Curie- Weiss 
G = — 244K, whereas Ref. [26| found J nn /&B = 200K. The former calculation agrees much better with experiment [27| 
which gave 8 — — 465K. Both groups studied the electronic band structure but it was not entirely clear what space 
group their calculations predicted. More detailed information on the symmetry of the structures comes from the 
first principles calculations of BF some of which included spin-orbit interactions. These calculations give a weak 
ferromagnetism of 0.18/zb per unit cell which is somewhat smaller than 0.3/xb per spin. [28] However, inclusion of spin- 
orbit interactions enabled BF to obtain the correct symmetry of the magnetoelectric behavior. Here we discuss in 
detail the symmetry properties of the various phases and experimental consequences such as the interactions between 
various order parameters (OP's) and the number and symmetry of the various domains which may be observed. 
Several of these issues were discussed by BF, but a more complete and systematic analysis is given here. 

Our approach to symmetry is similar to that of Ref. [2^ in connection with the Aurivillius compound SrBi2Ta20g 
(SBTO): we adopt the high-temperature tetragonal structure as the "reference" structure and analyze the instabilities 
at zero temperature which lead to the lower temperature phases. Although CMO and CTO do differ from SBTO, 
their crystal symmetry is the same as SBTO and hence many of the results we find here are similar to those for SBTO. 
Here we emphasize some of the experimental consequences of the symmetries we find (such as the enumeration of 
the different possible structurally ordered domains) and also we explore the nature of the macroscopic strains, the 
ferroelectric polarization, the magnetic ordering, and the coupling between structural distortions and these degrees of 
freedom. 

Briefly, this paper is organized as follows. In Sec. II we outline the basic approach used to analyze the symmetry 
of the systems in question. In Sec. Ill we give a symmetry analysis of the resulting phases which result from the 
structural instabilities found by BF. This analysis closely parallels that of Perez-Mato et al. 20]. In Sec. IV we discuss 
the second structural transition in which the other two irreps condense to reach the Cmc2i phase. In Sec. V we 
discuss the symmetry of the magnetic ordering. We also explore the coupling betweeen the distortions and the strains, 
the polarization, and the magnetic ordering. Throughout the paper we point out experiments which are needed in 
order to remove crucial gaps in our understanding of the structural phase diagram of these systems. In Sec. VI we 
enumerate the possible domains that can occur and briefly discuss their dynamics. Our conclusions are summarized 
in Sec. VII. 



II. SYMMETRY ANALYSIS 



Our symmetry analysis will be performed relative to the tetragonal I4/mmm structure which is stable at high 
temperatures. In this high-temperature reference structure one has atoms at their equilibrium positions R^°^(n, r), 
where 

R^{ni,n 2 ,n 3 ;T) = nxR lta + n 2 R2, a + n 3 R 3 , a + T a , 

where the lattice vectors are 

R x = (-1/2,1/2,1/2), R 2 = (1/2, -1/2, 1/2), 
R 3 = (1/2,1/2,-1/2) , 

where a labels Cartesian components, real space coordinates are expressed as fractions of lattice constants (so that, 
for instance, Ri denotes (— a/2, a/2, c/2)), and the number r labels sites within the unit cell at the position vector 
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r, as listed in Table HI The associated tetragonal reciprocal lattice vectors are 

Gi = (0,1,1), G a = (1,0,1), G 3 = (1,1,0), 
in reciprocal lattice units so that, for instance, Gi denotes 2n(0, 1/a, 1/c). 
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FIG. 2: (Color online) The group-subgroup structure arising from the first principles calculations. The transitions from I4/mmm 
are labeled with irrep that is condensing. In the final transition to Cmc2i we show below that the two remaining irreps condense 
to that all three irreps have condensed to form Cmc2i. This diagram does not deal with magnetic ordering. 



This work was stimulated by the first principles calculations of BF on the instabilities at zero temperature of the 
reference tetragonal system. For CMO the instabilities at zero temperature occur for the irreps X$ and X^ at the 
zone boundary X points which are K = (1/2, 1/2, 0) = qi and K = (1/2, —1/2, 0) = q 2 . (The superscript indicates 
the parity under inversion about the origin.) Since — q„ is equal (modulo a reciprocal lattice vector) to q n , the vectors 
qi and q2 exhaust the star of X. Also a zone center phonon of the two dimensional irrep T§ is nearly unstable. These 
results suggest the group-subgroup structure shown in Fig. [3J which is similar to that for SBTO420] We therefore 
consider structures having distorted positions given by 

R a (n,r) = R^(n,r)+u a (r)e lcl -^ Rknk) 

where q is the wave vector of the distortion mode, u(r) is the distortion taken from Tabled and e a p is the macroscopic 
strain tensor. The leading terms in the Landau expansion of the free energy of the distorted structure relative to 
the reference tetragonal structure will have terms quadratic in the strains e a p and the microscopic displacements 
u(t) within the unit cell. The fact that the free energy has to be invariant under all the symmetry operations of the 
"vacuum" (i. e. the reference tetragonal structure) restricts the microscopic displacements to be a linear combination 
of the basis vectors of the irrep in question. The basis functions are listed in Table Q] and the representation matrices 
for the generators of the irreps are given in Table HH29f in terms of the Pauli matrices 
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Using the basis functions w\_ and W 2 listed in Table Q] one can check that the matrices of Table [TTI do form a 
representation, so that for an operator O we have 



O 



Mu(O) Mia(O) 
M 21 (0) M 22 (0) 
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where the 'f's are normalized: 
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TABLE I: Basis functions *J Z) (R, r) = *^( r ) cos (q ' R ) and *2 Z) ( R > r ) = ^f,l( T ) ^(l ' R ) for the distortion vector under 
irrep Z, where Z = 5, 3, 2 indicates irrep F,T, , and X^ 1 ", respectively, and n = x,y, z labels the components. The <I/'s are 
normalized according to Eq. {T}. For the X irreps the wave vector of \pj^(R, r) is q^. For the wave vector is q = 0. For 
each site r we give the three components (n = 1, 2, 3) of the vector displacement. The values of the displacements in this table 
are not restricted by symmetry. The values of the structure parameters for CMO, taken from Ref. |7J, are p = 0.311, £ = 0.100, 
X = 0.205, and r = 0.087. 
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TABLE II: Representation matrices M^ 2 \ and M'" for the generators of the irreps -XT , Xp~, and , respectively. Here 

r' = Or and the cr's are the Pauli matrices. These matrices are related to those of Ref. 18'i 3Q| by a unitary transformation 
which, for takes (cr x , cr y , a z ) into (— er x , — cr y , cr z ). For X^ and X^ the unitary transformation takes (cr x ,(Ty,cr z ) into 
(— tr z , <T y , — cr x ). Unlike Ref. 18. we choose a representation in which the matrices representing the translations T„ are 
diagonal. 
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We now introduce OP's Q\ and Q 2 as amplitudes of these normalized distortions, so that a distortion <&z of the 
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irrep Z can be written as $z = Q\ '^f\ 
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and 
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We now interpret this as defining how the order parameters transform when the ^/n are regarded as fixed. Thus 
we have 
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Note that the OP's transform according to the transpose of the irrep matrices 
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III. UNSTABLE IRREPS 



A. X 3 - 

The irrep of the little group (of the wave vector) is one dimensional. However, since there are two wave vectors 
in the star of X, we will follow Ref. and construct the two dimensional irrep which incorporates both wave vectors 
in the star of X. The resulting two dimensional matrices are given in Table |TTJ From the basis functions for irrep 
X^ given in Table U one sees that the distortion can describe the alternating tilting of the oxygen octahedra about 
a (1, 1, 0) direction if q = qi and about a (1, 1, 0) direction if q — q2, as shown in Fig. [31 We now construct the form 
of the free energy when the distortion is given bv(3l| 

$ = Q 3 -( qi )vl4 3) +Q3-(q 2 )M4 3) . 




x = -10 1 2 x = -10 1 2 

q = (-1/2,1/2,0) = q 2 q = (1/2,1/2,0) = q % 

FIG. 3: (Color online) Schematic diagram of the displacements ^ 3 '(R, r), left and (R, t), right of the oxygen octahedra. 
All coordinates are in the parent tetragonal system. The arrows represent displacements within the x-y plane and the + and 
— signs, displacements collinear with the z-axis. The arrows represent the displacement of the apical oxygen at (0,0, x), which 
is the same as that at (0, 0, — \). The displacement of the shared apical oxygen at (0, 0, 0) is not shown. The algebraic sign of 
the displacements in the z direction of the equatorial oxygen at (x, y, r) in the unit cell is indicated. The displacement of the 
oxygen at (x,y, — r) is the negative of that at (x,y,r). The orthorhombic axes for Cmcm are X, Y, and Z. The distortion of 
left panel is obtained from that of the right panel by a 90° rotation about the positive z-axis. 

To construct the form of the free energy for this structure, note that wave vector conservation requires that the 
free energy be a function of Qa(qi) 2 and Qs(q 2 ) 2 because 2qi and 2q 2 are reciprocal lattice vectors but qi + q 2 is 
not a reciprocal lattice vector. Then, using the irrep matrices given in Table UH one can check that the free energy in 
terms of the X$ OP's must be of the form 

7YX 3 ) = ^(T-T 3 )[Q 3 ( qi ) 2 + Q3(q 2 ) 2 ] 

+ ^[Q 3 (qi) 2 + g 3 (q 2 ) 2 ] 2 

+«[Q3(qi)Q 3 (q 2 )] 2 + 0(Q 6 ) , (3) 
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where a > and T3 is the temperature at which this irrep becomes active if it is the only relevant OP). The first 
principles calculations of BF indicate that u > 0. To treat the coupling between the Q's an d the strains we consider 
the strain-dependent contribution to the free energy Fq e . The term whereby the Q's induce a strain is linear in the 
strain and thus we write FQ e as 

-^Qe = g Cnm£n€m + JklmtklQs (q?n) 2 

nm klm 
1 v 

— ^ / j ^nm^n^m ~\~ j 
nm 

where the first term is in Voigt[32[ notation, where 1 = x,x, etc. and only en, C12, C13, C33, C44, and cqq are nonzero 
under tetragonal symmetry and 

Vq, £ = ae X y[Qz (cu f - Q3 (q 2 ) 2 } + {P(e X x + e yv ) 

+7^][Q3(qi) 2 + Q3(q 2 ) 2 ] , (4) 

where a, /3, and 7 (and similarly below) are arbitrary constants. The sign of the shear deformation e xy depends on 
the sign of Q3~(qi) 2 — Q^i^) 2 , a result similar to that given in Ref. 

As the temperature is reduced through the value T3, a distortion of symmetry appears. For v > either one, 
but not both, of Qs(qi) and <5s(q2) are nonzero, so that e xy ^ and a detailed analysis shows that the resulting 
structure is Cmcm (#63). The first principles calculations of BF imply that v > 0, so, as indicated by Fig. [2j 
this possibility is the one realized for the RP systems we have studied. If v had been negative, then we would have 
had a "double-q" state (i. e. a state simultaneously having two wave vectors) with |Q3(qi)| = |Qa(q2)|, so that 
e xy = and the distorted structure would be the tetragonal space group P4 2 /mnm (#136). These two space groups 
are listed in Ref. [18| as possible subgroups of I4/mmm which can arise out of the irrep -^3". These directions in 
OP-space [Qs(qi) — Q3(q2)-space] are stable with respect to perturbations due to higher order terms in Eq. 
which are anisotropic in OP-space as long as T is close enough to T3 so that these perturbations are sufficiently 
small. We note that Ref. [TH lists Pnnm (#58) as an additional possible subgroup. This subgroup would arise if v 
were exactly zero, in which case, even for T arbitrarily close to T3, the OP's would be determined by higher order 
terms in Eq. ([3]) and would then not be restricted to lie along a high symmetry direction in OP space. However, to 
realize this possibility requires the accidental vanishing of the coefficient v and the analogous sixth order anisotropy. 
This possibility should be rejected unless additional control parameters, such as the pressure or an electric field, are 
introduced which allow access to such a multicritical point. 10] As of this writing, the Cmcm phase (or any other phase 
intermediate between 14/mmm and Cmc2i) has not been observed, although it has been observed in the isostructural 
systems LaCa2Mn 2 07[22| and Bin aaC&i ^KMn9Q7|23l| 

For future reference, it is convenient to introduce orthorhombic axes, so we arbitrary define them as in Fig. [4] If 
we condense Q3(qi), the tilting is about Y, the tetragonal [1 1 0] direction, whereas if we condense Q3(q2), the tilting 
is about Z, the tetragonal [1 1 0] direction. The sign of the OP indicates the sign of the tilting angle. Note that 
there are four domains of ordering, depending on which wave vector has condensed and the signs of the OP. These 
correspond to the four equivalent initial orientations of the tetragonal sample. In M3V2O8, such a domain structure 
for two coexisting OP's was confirmed experimentally. pjij and it would be nice to do the same here. In Sec. VI we 
give a more detailed discussions of domains. Here we note that changing the sign of the OP Q3(qi) is equivalent to 
a unit translation along x or y. Thus, as for an antiferromagnet, the phase of the order parameter within a single 
domain has no macroscopic consequences. 

Since the Cmcm structure has a center of inversion symmetry, it can not have a nonzero spontaneous polarization 
and an interaction Vq.p which is linear in the polarization does not arise. It may seem mysterious that when we 
introduce a distortion which is odd under inversion, we still have a structure which is even under inversion. The point 
is that the distortion due to X^ is odd under inversion about the origin, but is even under inversion with respect 
to (1/2,0,0), about which point the parent tetragonal structure also has a center of inversion symmetry. To see 
the inversion symmetry of the distortion about (1/2,0,0) from Fig. [3l note that inversion displaces the arrows and 
reverses their direction. Since a plus sign represents a positive z distortion at z — t and a negative z distortion at 
z = — t, inversion takes a plus into a plus and a minus into a minus. 



As in the case of X 3 we construct the two dimensional irrep which incorporates both wave vectors in the star of 
X, whose matrices are given in Table |Hj The allowed basis functions are given in Table U and are represented in 
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FIG. 4: (Color online) The tetragonal axes (lower case) and the axes (capitals) we use to describe the orthorhombic phases. 
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Fig. [5j Once the configuration for one wave vector is determined, that of the other wave vector follows from a 90° 
rotation about the z-axis. So these are two different configurations (one for each X wave vector) which have the 
same free energy. These two configurations differ in their stacking (which takes point O into point O'). This stacking 
degeneracy means that this transition takes the I4/mmm structure into one of two different settings of the space 
group which we identify below as Cmca. One sees that in this distortion the oxygen octahedra are rotated about 
the crystal c axis as if they were interlocking gears. Notice that in the left panel of Fig. [5] all the clockwise turning 
octahedra have r x > r y and all the counterclockwise ones have r x < r y , whereas in the right panel the clockwise 
turning octahedra have r y > r x and the counterclockwise ones have r y < r x , where r x is the radius of the elliptically 
distorted octahedra along x (or nearly x) and r y is the radius along y. Symmetry does not fix the sign of the radial 
distortion. Changing the sign of the radial distortion would lead to a different (inequivalent) structure in which the 
sign of all the radial arrows for both wave vectors would be changed. It would be interesting to observe this radial 
distortion in either (or both) calculations and experiment. 

As before, we introduce OP's Q2(qi) and Q2(q2) by considering the free energy for which the distortion from 
tetragonal is given by Q2(qi) times the distortion for qi plus Q2(q2) times the distortion for q 2 . As in Eq. @, 
the transformation properties of these OP's are determined by the matrices of Table [TTJ As before, wave vector 
conservation requires that the free energy be a function of Q2(qi) 2 and Q2(q2) 2 - Using the transformation properties 
of the OP's, we find that 



where a' > 0, u' > 0, and T2 is the temperature at which this irrep would become active, if it were the only relevant 
irrep. The free energy of the coupling to strains is 



where a', /3', and 7' are arbitrary coefficients. The coupling to the shear strain is allowed because both e yx and 
Q2(qi) 2 — Q2(q2) 2 are odd under 7?4 and even under nid- The term [e xx — e yy ] [<22(qi) 2 — Q 2 (q2) 2 ] is not allowed 
because it is not invariant under m c i- Since this structure is even under inversion it can not have a spontaneous 
polarization, so -Fq,p, the interaction linear in the polarization is zero. 

Now we discuss the phases which Landau theory predicts. If v' > 0, then either Q2(qi) = or Q2(q 2 ) = 0, so that 
€ xy 7^ and a detailed analysis shows that the resulting structure is Cmca (#64). Similarly, if v' < 0, then one has 
a double-q state with |Q2(qi)| = |Q2(q2)|, the orthorhombic distortion vanishes, and we have the structure P4/mbm 
(#127). The other possibility listed in Ref. [l8|is Pbnm (#35). As before, we ignore this possibility because it requires 
the accidental vanishing of v' . The first principles calculations of BF indicate that v' > as required by Fig. [21 



J"(X 2 ) = ^(T-T 2 )[Q 2 (qi) 2 + Q2(q2) 2 ] 

+ i7i'[Q 2 (qi) 2 + 02(q 2 ) 2 ] 2 
+«'[Q 2 (qi)Q 2 (q2)] 2 +0(Q 6 ) , 



= a 



*W<52(qi) 2 - Q2(q2) 2 ] + \P' (ex* + e vy ) 
+7'^][Q 2 (qi) 2 + Q2(q 2 ) 2 ] , 
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q 2 = (-1/2,1/2,0) q 1 = (1/2,1/2,0) 

FIG. 5: (Color online) Schematic diagram of the displacements for X% , which are confined to the x-y plane. All coordinates 
are in the parent tetragonal system. The large black arrows give the rotational distortion. The much smaller anisotropic radial 
distortion (which is allowed in X£) is shown in magenta. In both panels the distortions are shown for sites in the planes z — ±r 
and z — 1/2 i t. The distortion is an even function of r. The distortion of the right panel is obtained from that of the left 
panel by a 90° rotation about the z-axis. The translation vectors in the x-y plane are (1,1,0) and (1,-1,0). In the left panel 
the translation vector Ar (which takes O into O' and obeys exp(iq • Ar) = 1) is (1/2, 1/2, 1/2) and in the right panel it is 
( — 1/2, 1/2, 1/2). These two wave vectors thus give rise to the two settings of the side-centered orthorhombic lattice. 



c. r 5 - 

The allowed distortion of the two-dimensional irrep is given by the basis functions of Table U and this distortion 
breaks inversion symmetry. As before, we assume that the distortion >F from I4/mmm is given by <5s,i times the 
distortion "J/j 5 plus Q^,2 times the distortion for ty 2 > where the 'S/'s are given in Table Q] and are seen to transform 
simlarly to x or y. It follows that the OP's Qs,i and (55,2, transform as x and y, respectively. So under tetragonal 
symmetry the expansion of the free energy F of the distortion 'F in powers of Qs,i and (55,2 assumes the form 

^(r 5 ) - y(T -T 5 )(Ql 1 + Ql 2 ) + ±u"(Ql 1 + Ql 2 ) 2 

+v"Ql^Ql 2 + 0{Q & ) , (5) 

where a" > 0, u" > 0, and T5 is the temperature at which irrep becomes active (if it were the only relevant irrep). 
Since BF find that is not actually unstable for CMO, T5 < 0. Using the fact (see Tables HI or HT|) that Q51 and 
(55,2 transform like x and y, we see that 

V Qt = ae zz {Q\ x + Ql y2 )+p{e xx + e vv ){Q% x + Q\ 2 ) 

+j(e X x ~ £yy)(Qs,i ~ Qla) + 5e xy Q b ^Q ba ■ (6) 

Because this phase is not centrosymmetric, it can support a nonzero spontaneous polarization. This will be discussed 
in a later subsection. 

The free energy of Eq. gives rise to two principal scenarios. If v" > 0, then either Q51 or Q5 2 (but not both) 
condense at T5 . Then Eq. JS]) indicates that e xy = and e xx 7^ e yy and a detailed analysis of the distortions indicates 
that we condense into space group Imm2 (#44). Alternatively, if v" < then Eq. ([5]) indicates that | 1 1 = IQ5.2 1 - 
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Equation (jSJ) implies that e xy 7^ and e xx = e yy and we condense into space group F2mm ($42). As before, these 
distorted space groups agree with the results in Ref. flSl . However, an additional subgroup is listed as realizable from 
this irrep, namely Cm (#8). As before, to realize this possibility requires the accidental vanishing of v", a possibility 
we reject. The first principles calculations of BF imply that v" < so that space group F2mm would be realized if 
were to condense first and both OP's would have equal magnitude. There are then four possible domains of ordering 
corresponding to independently choosing the signs of Qs,i an d Q5,2- 

IV. COMBINING IRREPS 

Now we consider what happens when we condense a second irrep. Although we assume that X$ is the first irrep to 
condense, our discussion could be framed more generally when the ordering of condensation of the irreps is arbitrary. 
We assume a quadratic free energy T-z of the form 

= §(T - 7 1 3 )[g 3 (qi) 2 + Q 3 (q 2 ) 2 ] 

+^(T-T 2 )[Q 2 (q 1 f + Q 2 (q 2 ) 2 } 

+^(T-T 5 )(Y 2 + Z 2 ) , 

where Y§ = X^i — Xs^, Z$ = X5.1+X5.2, we assume that T3 > T 2 > T5. For the RP systems wave vector conservation 
and inversion invariance implies that the cubic terms in the free energy must be of the form 

V c = Q2(qi)Q3(qi)[rQ B ,i+*Q5,2] 

+Q2(q. 2 )Q 3 {c[ 2 )[sQ 5A + uQ 5a ] . 

To make this invariant under rrid we require r = t and s = —u. Invariance under 72-4 leads to s — r, so that we may 
write 

V c = -r[Q 2 (q 1 )Q 3 (q 1 )Z 5 +Q 2 ( q2 )Q 3 ( q2 )Y 5 } . (7) 

As the temperature is reduced, Q 3 is the first OP to condense, with either Qa(qi) 7^ or Q 3 (q2) 7^ 0, as dictated by 
the quartic terms we considered previously in Eq. ([3]). Then, assuming Q 3 (qi) has condensed, we have effectively 

Fa = ™(Q 3 (qi))Q 2 (qi)Z 5 = w'QMZs , 

where (<5 3 (qi)) indicates the value of Q 3 (qi) which minimizes the free energy. Thus the Q 2 and Q5 variables are 
governed by the quadratic free energy 

^Q 2 ,z = i(T-T 2 )Q 2 ( qi ) 2 + i(T-T 5 )Z 5 2 

+w(Q 3 (q 1 ))Q 2 (q 1 )Z 5 + \(T- T 2 )Q 2 (q 2 ) 2 

+i(T - T 5 )Y 5 2 + w(Q 3 (q 2 ))Q 2 (q 2 )Y 5 . 

Suppose that Q 3 (qi) is the first variable to condense. The effect of the cubic term is to couple <52(qi) and Z$ so 
that the variable that next condenses as the temperature is lowered is a linear combination of Q 2 (qi) and Z5. If w' 
is small compared to T 2 — T5, then the new transition temperature will be approximately T 2 w T 2 + w' /(T 2 — T5) 
and the condensing variable Q2(qi) will dominantly be Q2(qi) with a small amount of Z$ admixed to it. The 
important conclusion is that we have two families of OP's: {Q 3 {qi),Q 2 {qi), Z 5 ) and {Q 3 (q 2 ),Q 2 (q 2 ),Y 5 ). At the 
highest transition one OP (we assume it to be Q 3 ) of one of the two families will condense. At a lower temperature 
the two other OP's of that family will condense. Independently of which OP first condenses one will reach one of the 
equivalent domains of the same final state, as Fig. [5] indicates. There are eight [Hf such equivalent domains because 
we can independently choose between a) the wave vectors qi and q 2 , b) the sign of Q 3 (q), and c) the sign of Q2(q)- 
These domains are discussed in detail in Sec. VI. For simplicity, we do not extend this analysis to include several 
copies of the various irreps involved. The symmetry of the displacements when the irreps of the qi family are present 
is shown in Fig. [5] 
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The cubic term of Eq. (19) guarantees that when Q 3 (qn) condenses, the lower-temperature transition always involves 
the condensation of Q 2 (qn) and the appropriate T§ OP. We selected this cubic interaction to be dominant in order 
to arrive finally at the observed Cmc2i phase, as discussed in Appendix [5] There is also a cubic interaction of the 
form 

V'c = r[Q2(m)Q 3 (q2)V(q M )+Q2(q2)Q3(m)W(q M )} , 

where q^f = (1,0,0) and Vm and Wm are operators that transform according to the appropriate two-dimensional 
irrep (M^~) so that V' c is an invariant. This interaction (in contrast to Vc) involves a doubling of the size of the unit 
cell at T2. Since no such doubling has been seen, we conclude that V c is dominated by Vc- 



A 




FIG. 6: (Color online) The symmetry of the displacements of the equatorial oxygen ions at tetragonal sites (±1/2,0, r) and 
(0, ±1/2, r) and that of apical oxygens (e. g. e) at (0, 0, x) in the Cmc2i phase. Arrows represent displacements in the xt-yt 
plane and the + and — signs those in the z t direction. Each octahedron (e. g. a, b, c, d, e) is surrounded by a dash-dot 
oval. The tetragonal (orthorhombic) unit cell is bounded by the dashed (full) rectangle. The spontaneous polarization P and 
the wave vector q are shown at the right along with the tetragonal and orthorhombic coordinate axes. The displacements of 
irrep T§ (which contribute to P ^ 0) are in parentheses. The other symbols have the following meaning. The octahedral 
rotation (the largest arrows) about zt with the smallest arrows for the radial distortion come from irrep X^ and that about 
Zo (the apical arrows, the ±'s, and the -'s) comes from irrep . The only symmetries remaining in the Cmc2i phase are 
m z = mx = ±1 and a glide plane whose mirror is the dashed line and whose displacement is the dashed arrow. The glide 
operation takes site 1 into 2 and 3 into 4. 



It is interesting to consider the mean-field temperature dependence of Q 3 (qi), for example, from the free energy 

T 



\{T -T 3 )Q 3 ( qi ) 2 + \uQMit 



For T slightly below 7a one has |Q 3 (qi)| 
find for T < T 2 that [36] 



+wQ 3 (q 1 )Q 2 (qi)[Q 5 ,i + Q5,2] ■ 
(T 3 - T) 1 / 2 . For T near T2 we treat the term in w perturbatively and 



Qs(qi) 



Slightly below T2 the variables Q2(qi) and (Qs,i 
to (f 2 - T) 1 / 2 . Thus the second term in Eq. (0), 
results are illustrated in Fig. [7] 



[(a/u)(T 3 - T)] 1/2 

wQ 2 (qi)(Q54 + Qs,2) 
2a(T 3 - T) 



(8) 



±Q5,2) are proportional (with different constants of proportionality) 
which is only nonzero for T < T2, is proportional to (T 2 — T). These 
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P,Q's 




FIG. 7: (Color online) Temperature dependence of Qz for T near T3 and (lower curve) of Qz, Q2, and Qs for T near T2, as 
discussed below Eq. ([8]) and of P as described by Eq. (flO|) . 



One can similarly analyze the temperature dependence of the strains. From Eq. @ we see that the strain e xy (in 
tetragonal coordinates) is zero for T > T3. Just below T3 one has 

e xy = a[Q3~(q 2 ) 2 - Q3 (qi) 2 ]/ 'c 4 4 ■ 

Since Q3 oc (T3 — T) 1 / 2 , this indicates that just below T3 one has 

e^y = ±( T 3 - T) , 

the sign depending on which Q3(q n ) has condensed. In addition, the coupling between strains and the orientational 
OP's indicate that at the structural transitions there will be a jump in slope of the diagonal strains e aa . 



V. COUPLING TO DIELECTRIC AND MAGNETIC ORDER 



In the following two subsections we consider the coupling between the lattice distortions and a) the spontaneous 
polarization and b) magnetic long range order. 



A. Dielectric Coupling 



The dielectric free energy is 



where Vq is the coupling with the distortion modes and which is linear in the spontaneous polarization P. This 
coupling will be zero for a structure which has a center of inversion symmetry. is the only irrep which, by itself, 
breaks inversion and for it we have 

This follows because Q5.1 transforms like x and Q5.2 transforms like y. Using this transformation property, we infer 
from Eq. ([7]) a contribution to Vq of the form 

^ 2 ' 3) = -A , [Q 2 (qi)Q a (qi)(i , x+i , v) 
+Q2(qa)Qs(qa)(i , x • 

These couplings indicate that 

P x = AxEQ5,i + A'x E [Q 2 (qi)Q3(qi) + Q2(q2)<33(q2)] 

Py = Ax£Q5.2 + A'xB[Q 2 (qi)<23(qi)-<32(q2)Q3(q2)] ■ 

(9) 



11 



We may simplify this by minimizing the free energy of Eqs. (JS|) and to write 



Qs,i = ^^r(Q2(qi)Q3(qi) + Q 2 (q 2 )Q3(q2) 



Q5,2 = _ _ ( Q2(qi)Q3(qi) - Q2(q2)Q3(q2) 



T-T 5 ^ 
Thus 

p x = r^Q 2 (qi)Q3(qi) + Q2(q 2 )g 3 (q2)^ 

Py = r^Q 2 (qi)Q 3 (qi)-Q 2 (q 2 )Q3(q 2 )^ , (10) 

where r = xe [A'+Ar/(T— T 5 )] . As shown in Appendix[B] the dielectric constant will have a small amplitude divergence 
at a temperature T = T2 near T2, as is often seen in systems with magnetization induced polarization. [37||38| . Note 
that there are two mechanisms for a spontaneous polarization proportional respectively to A' and A. The term in 
A may be viewed as being the polarization due to displacement of the charged ions in the polar irrep . This 
displacement is induced by the presence of the other two OP's Q 2 and Q3. The term in A' is due to the modification 
in the electronic structure proportional to Q2(q)Q3(q) when the displacement due to is zero. The numerical work 
of BF indicates that there is a large (~ 5/zC/cm 2 ) spontaneous polarization which arises when Q5 is zero and therefore 
indicates that mechanism b) is dominant. 

For T slightly below T2, Q3 is essentially constant and the other Q's are proportional to (T2 — T) 1 / 2 within mean 
field theory, as illustrated in Fig. [Jj Equation ([TO]) indicates that P is parallel to q and is proportional to Q2Q3, so 
that near T2 one has the mean- field result that |P| ~ (T2 — T) 1 / 2 , as shown in Fig. [7J From the work of Mostovov. [391 
magnetically induced polarization was not expected to have P || q, although P || q, was found experimentallvpof 
and explained from general symmetry arguments. |40l)4ll| Equation (|9]) indicates that a polarization will be induced 
parallel to q along one of the four (1, 1, 0) directions according to the signs of the OP's as indicated by Eq. (j9|). 



B. Magnetism 



We now discuss the magnetic structures that can appear in this system. Ca3Mn2 07 becomes antiferromagnetic 
at Tjy — 115K. 14] For this discussion we will work relative to the parent tetragonal lattice but will introduce the 
orthorhombic coordinates for the spin vectors so that 

(Sx)o = (S z ) t , (S Y )o = (S x ) t + (Sy) t , 
(Sz)o = (S x )t — {S y )t ■ 

In the simplest approximation the antiferromagnetic structure of a single biliayer (consisting of one layer of Mn ions 
at z t = x Q = £ and another at z t = x Q = — £) is that of two square lattice antiferromagnets stacked directly on 
top of one another [sites at (0,0, £) and (0,0, — £)], so that all near neighbor interactions proceed via nearly 180° 
antiferromagnetic Mn-O-Mn bonds. The first principles calculations of BF and the data from Ref. [3 indicate that 
dominantly the spins are perpendicular to the plane. This structure is shown in Fig. [8l Therefore we assume that the 
dominant order parameter is G(q), where we use the Wollan-Koehler[42j symbols to represent the configurations of 
a single bilayer shown in Fig. [9j In the G z = Gx configuration one can have either q = qi or q = q2. The difference 
between Gx(qi) and Gx{Q2) is m the different phase of one bilayer relative to the other. In any case all 5 magnetic 
nearest neighbors of a central spin are oriented antiparallel to it. 

We will accommodate the following structures based on the parent tetragonal lattice. We have choices for the wave 
vector, namely q = 0, q = qi and q2 and S a (0, 0, £) = ±S a (0, 0, — £). So for each component of spin we have six 
candidate structures. If q = and S a (0, 0, — £) = 5(0, 0, £), then we have the "F" (ferromagnetic) structure. If q = 
and 5 Q (0,0, — £) = —5(0,0, £), then we have the "A" structure shown in Fig. |H1 If q = q n , then each plane of the 
bilayer consists of a square lattice antiferromagnet. The two planes of the bilayer can be coupled either so that adjacent 
spins in the planes are parallel (this is the "C" structure) or so that they are antiparallel (this is the "G" structure), 
as shown in Fig. |H1 The C and G structures each come in two versions depending on whether q = qi = (1/2, 1/2, 0) 
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FIG. 8: (Color online) A single bilayer of Mn ions when the spins are assumed to lie in the plane. This is a Gx configuration. 
The stacking of adjacent bilayers (not shown) is such that they are displaced transversely by (1/2,1/2,0) relative to one 
another. The phase of antiferromagnetic order in adjacent bilayers is determined by the wave vector which is either qi or q2. 
This configuration of pseudovectors is odd under inversion about the origin (x = y — z — 0). 




FIG. 9: The A, C, F, and G configurations of spin components for a single bilayer in the Wollan-Koehler [42l] scheme. 




A C 



FIG. 10: Spin states of a bilayer of Mn ions. The plus and minus signs represent the signs of any component of spin. Left: 
the "A" configuration. Right: the "C" configuration. The circled symbols represent the spins in planes at z — 1/2 + £ and 
z = 1/2 — £ for q = qi. For q = q2 the circled + and circled — signs of the C configuration are interchanged. 



or q = q2 = (1/2, —1/2,0), as discussed in the caption to Fig. [TUJ We now introduce orthorhombic coordinates, so 
that Sx = S z , Sy — S x — S y , and Sz — S x + S y . Then we have the symmetries given in Table Hill 
The magnetic free energy is F = Fq + V, where Fq is the free energy of the G structure: 



F a = 



(T -T N + K a )Y,G a {cinY 
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TABLE III: The symmetry for components X, Y, and Z of the spin, a pseudovector. Here Sx = S x , Sy = S x + S y , Sz = S x ~S y , 
where capitals refer to orthorhombic and lower case to tetragonal. 
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qi 


+ + + 




- + 


+ 




C 


q2 


+ + + 




- + 


+ 





-a'^G ct (q„) 2 Q 3 (q„) 2 



0[G(q„) 4 



and 



V = 



[i a Fl + ^C Q (qi) 2 +r a A 2 c 



(11) 



The quartic term in Fq is such as to ensure that the magnetic ordering vector is the same as that of the paramagnetic 
structure. Now we want to see what other magnetic OP's are induced by the condensation of the dominant G ordering 
in the presence of the nontetragonal distortions. The magnetoelastic interaction we invoke has to be quadratic in the 
magnetic variables in order to be time- reversal invariant. So we consider a cubic potential which contains terms of 
the form 

Gx{q.n)K a Qp , 

where K is C(q n ), F, or A and Q is (q n ), Q$ (q„), Z 5 = Qs,i +<?5,2, or Y 5 = Q 5i i - Q 5 , 2 - We start by considering 
only terms involving qi. The terms involving q2 will later be obtained from those involving qi by applying the 
four- fold rotation IZ4. The only terms which are consistent with inversion symmetry and wave vector conservation 
are those of the form 



G x (q 1 )A a Q+(q 1 ) , G x ( qi )F a Qs (<n) , 
G x (( ll )C a (q 1 )Z 5 . 

(Y 5 = for q = qi.) We now use Table Hill to require invariance under and m z , so that the interaction which has 
the correct symmetry is 



aG x (m)F Y Q^ ( qi ) + 6G?x(qi)C z (qi)Z, 



5 • 



(12) 



Now we use 



TliGx (q«) 

K.4Z5 



Gx{<l3-n) , TZ^Fy — —Fz , 

-CV(q 2 ) , (qi) = Q^"(q 2 ) 



Thus, in all, the lowest order magnetoelastic coupling Vmq is 



V MQ = aGx(qi)FyQ 3 -( qi )+6Gx(qi)Cz( qi )Z 5 

-aGx^FzQs (q 2 ) - 6Gx(q 2 )Gy(q 2 )y 5 ■ 



(13) 



To see what this means, it is helpful to recall that qi (q 2 ) lies along the orthorhombic Z (Y) direction. Thus the weak 
ferromagnetic moment F is perpendicular to q. Note that the wave vector is already selected as soon as tetragonal 
symmetry is broken. Say qi is selected. Then, in addition, the sign of Q3 (qi) was also selected when tetragonal 
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symmetry was broken. Then, when magnetic long-range order appears, it can have either sign of Gx(qi), but the 
sign of Gx(qi)Fz is fixed by the interactions within the system. 

Assuming the G configuration to be dominant and using Eq. (jlip . we thus have two scenarios. If q = qi, then we 
have the magnetic OP's 



[G x (qi),F Y ,Cz(<li)] 
F Y _ oQaCqi) 



Gx(qi) 

If q = q2, then we have the magnetic OP's 

[Gx(q 2 ),C r (q 2 ,F z )] 
Gy(q 2 ) = W5 
G x {q.2) v z 



with 

Cz(qi) 

Gx(qi) 



(14) 



with 



G x (q 2 ) 



gQ 3 (qg) 

Mr 



(15) 



In all the above results, since the cubic coupling combines Z$ with Q^(qi)Q^(qi) and similarly for q = q2, we should 
replace Z§ by a linear combination of Z$ and Q3 (qi)Q^" (qi) and Y5 by a linear combination of Y5 and Q3 {^xjQ^isti)- 
The result of Eqs. (fT4")l and (fT5")) agrees with the magnetic structure determination of Ref. [H] and with the symmetry 
analysis of BF, except that here we emphasize the relation of the magnetic ordering to the preestablished wave vector. 
These results are summarized by Fig. [TT1 





FIG. 11: (Color online) The equilibrium orientation of the staggered magnetization (N), the polarization (P), and the weak 
ferromagnetic moment (F) for q = qi (left panel) and for q = q2 (right panel). The axes along which these vectors lie are 
fixed by symmetry. However, one can find domains in which any (or all) these vectors can be reversed. For instance, P is 
proportional to Q2Q3 whereas F is proportional to Q3. 



It is interesting to note the possibility of switching the direction of the polarization (magnetization) by application 
of a sufficiently strong magnetic (electric) field. For this discussion it is useful ro refer to Fig. [11] Suppose the sample 
initially has condensed wave vector qi. (See the left panel of Fig. [TTJ) Applying a large enough magnetic field in the 
Z direction (i. e. parallel to qi) will cause a reorientation of the magnetization of the system so that Fy = and 
Fz =/= 0. Thus Vmq indicates that the system will make a transition from wave vector qi to wave vector q2. Since the 
direction of the polarization is tied to the direction of the wave vector, this transition will also involve a rotation of 
the polarization. Likewise, suppose the sample initially has the polarization and the condensed wave vector collinear 
and parallel to qi . Equation (fl4l indicates that the ferromagnetic moment is along Y (which is perpendicular to qi ) . 
Then, if one applies a large enough electric field parallel to q2, then the polarization will rotate into the direction 
of q2 and since the polarization and wave vector are constrained to be collinear, the wave vector will now be along 
q2- Then Eq. ([T5]) indicates that the weak ferromagnetic moment will be rotated into the Z direction (which is 
perpendicular to q2). However, the applied fields required to accomplish these switchings may be extremely large in 
view of the structural reorganization involved. 

Since F a transforms like a pseudovector, one has symmetry-allowed interactions with F a replaced by H a where H 
is the applied magnetic field. So from Eq. (fT5)) we have a magnetic-field dependent contribution to the free energy of 
the form 

Vk - ciGjc(qi)Qi"(qa)fly + c 2 Gx(q 2 )Q^ (q 2 )#z , (16) 

which indicates that when tetragonal symmetry is broken (so that Q3 is nonzero), the magnetic field acts like a 
field conjugate to the antiferromagnetic OP Gx- Consequently, dGx/dH a will diverge as the lower transition is 
approached for a — Y or a — Z, according to which wave vector has condensed. This suggests a neutron scattering 
experiment to measure Gx near the lower transition as a function of the magnitude and direction of H. 
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VI. DOMAINS 



Here we discuss in more detail the possible domain structures and give a brief discussion of the dynamics of domain 
wall motion. We first enumerate the various domains of order parameters which can exist as the temperature is 
lowered through the various phase transitions. As a preliminary one should note that within a single domain the 
phase of an OP at wave vector qk can not be experimentally determined. However, if more than one such OP is 
present, then their relative phases can be accessed experimentally. In the discussion that follows we will determine 
the phase of the OP's relative to that of Q3 which is not determined. As the temperature is lowered through T3, four 
possible domains are created, with the choices of sign of Q3 and the two choices for the wave vector, as shown in 
Fig. [TJ] The value of the wave vector within a single domain is experimentally accessible via a scattering experiment. 
Although the phase of Qs(q) can not be established within a single domain, the fact that different such domains 
do exist can be established by observation of a domain wall separating domains having the same values of the wave 
vector. Such an experiment to observe a so-called phase domain wall has recently been done in another system. (43j 
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FIG. 12: Domains in CMO. At T3 the value of the wave vector q and of Qa(q) are selected. At each succeeding transition a 
two-state order parameter condenses to further break symmetry. The domains for positive Qa(q) are macroscopically indistin- 
guishable from those for negative Qa(q) because they differ by a unit translation in the x-y plane. The signs of the OP's P and 
the weak ferromagnetic moment M depend on the constants r of Eq. (ffUJ and a of Eq. (I12|l . For CTO (which is nonmagnetic) 
the section at and below TV does not apply. 



Next, as the temperature is lowered through the lower structural transition at T = T2, the OP Q2 is condensed 
(and, as we have seen, is accompanied by Q$), giving a total of eight domains, four with Q 3 > and four with 
Q3 < 0. These sets of four differ from one another only in the inaccessible phase of the OP Q3. However, the four 
domains having a given sign of Q3 can be distinguished from one another, since they correspond to the two choices 
of wave vector (which is easily experimentally accessible) and the two choices of sign of Q2 which leads to distinct 
orientations of the spontaneous polarization (which is also easily experimentally accessible). 

Finally, when the temperature Tn is reached, the OP Gx(q) which describes the antiferromagnetic order is also 
accessible because it is coupled to the weak ferromagnetic moment. In that way we can identify the 16 domains up to 
an uncertainty in the sign of Q3. As we indicated above, the uncertainty is, in principle, accessible in that the phase 
domain wall can be observed. Note that the transformation (Q3, Q2, G x ) — > (— Q3, — Q2, —Gx) leaves the observables 
P, M, and q invariant. 

The above discussion assumes the existence of domains which, unlike ferromagnetic domains, do not have an obvious 
energetic reason to exist for T3 > T > T2, where P is zero. In particular, the phase domains, if they exist, would 
appear below the upper ordering temperature T3, where Qs(q) condenses, in which one has, for a given wave vector 
(which can be selected by applying a suitable shear stress), the two possible choices of sign of Qs(q). 
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FIG. 13: Displacements of the equatorial oxygens as in Fig. [3] An antiphase domain wall in the phase of a tilted octahedral 
phase in an n — 1 RP system. The phases to the left and to the right of the domain wall are indistinguishable since they only 
differ in the phase of the ordering. In this regard this situation is analogous to domains in an antiferromagnet. If the domain 
wall were a perfect plane, there would be no distortion energy of the octahedra. However, when the wall has a kink, as shown, 
then there are two possibilities. First, the tilting a, b, c, d can heal to zero, but this distortion will involve an intraoctahedral 
distortion of high energy. Secondly, it is possible that a, b, c, d (or analogously e, f, g, h) will alternate in sign, end at an 
antikink, so that no intraoctahedral energy is involved, but the tilting of the octahedra is not the lowest energy tilting of the 
established phase and the energy of the kink is proportional to the distance between the kink and the antikink. 




FIG. 14: (Color online) Domain walls in CMO for Tz > T > T%. Here one has low energy 90° walls between phases A and B, 
B and D, and C and D. Antiphase regions A and D are separated by two such walls. If the distance between regions A and D 
is not too large, it may be possible to establish the existence of phase domains using the technique of Ref. |43j. 



Upon cooling through T 3 it seems likely that one would have different nucleation sites from which ordering would 
develop. As a result, one would have randomly chosen signs of the OP Qs(q) in different regions of the sample. The 
question then arises, would these domains coarsen and the sample then become a single domain sample? To answer 
this question we need to study the energetics of domain walls. This involves understanding structural defects, as were 
studied in Refs. and For this discussion, I first consider the simpler case of an n = 1 RP system for which 
the tilt is around a [100] direction. In Fig. Q2] a domain wall is shown separating two phases to the left and right of 
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the wall which differ only in their phase. One sees that if the domain wall is a perfect plane, then there will be no 
intraoctahedral energy involved in the wall. Indeed, the only energy of the domain wall will be that from the potential 
that keeps the plane in its location. However, as explained in the caption to Fig. 1131 kink formation requires a large 
energy. Notice the difference between domain walls in this system and those in an Ising antiferromagnet. In an Ising 
antiferromagnet the domain wall energy is proportional to the length of the wall. Having a kink simply increases the 
energy by one unit of exchange energy. Here a kink either has the nonlocal string-like energy of a kink-antikink or its 
has the large energy needed to deform the octahedra. Thus, it seems that domain wall motion may be inhibited and 
therefore it is possible that phase domains, once created will remain in the sample. 

From these examples one concludes that a domain wall tends to form in planes of "minimum contact," i. e. in 
planes which intersect the least number of shared oxygen ions. This is why [100] walls are preferred over [110] walls 
for systems, such as CMO, in which tilting occurs about a [110] direction. [44j Of course, [001] walls (stacking faults) 
probably have the least energy. 

VII. CONCLUSIONS 

In this paper we have explored the rich structure of structural, magnetic, and dielectric ordering in the Ruddlesden- 
Popper compound CaaM^O?, using Landau theory to analyze symmetry properties. Our approach is similar to that 
used by Perez-Mato et al. |20| to study the Aurivillius compounds. Most of the symmetry relations we find are 
explicitly corroborated by the first principles calculations of Benedek and Fennie. An important aspect of our work 
is to motivate a large number of experiments which can elucidate the relations between the various order parameters. 
Specifically, we summarize the conclusions from our work as follows 

• 1 The most important aspect of our work is that we introduce order parameters (OP's) for all the irreducible 
representations (irreps) for all the wave vectors of the star which is active in the ordering transitions. The OP's 
describe distortions from the parent high symmetry tetragonal lattice which exists at high temperature. This enables 
us to discuss the induced (nonprimary) OP's such as the spontaneous polarization, the weak ferromagnetism, and the 
clastic strains. 

•2 In conformity with established results [l8| (but rejecting multicritical points) we treat the group-subgroup struc- 
ture obtained by the first principles calculations of Benedek and Fennie. [l9| We give an OP description of the sequence 
of structural transitions which is predicted for CaaM^Oj and Ca3Ti207, namely I4/mmm — > Cmcm — > Cmc2i the 
same as that found for a similar perovskite by Perez-Mato et aL|20| 

• 3 The ordering involves two families of domains, one family for each of the two X wave vectors. At the lower 
structural transition, a ferroelectric polarization appears parallel to the wave vector. Below that there is an inde- 
pendent magnetic ordering transition to an antiferromagnetic state in which the stacking of the magnetic bilayers 
depends on the wave vector which was selected when the tetragonal symmetry was broken. A weak ferromagnetic 
moment develops perpendicular to both the staggered magnetization and the polarization. 

•4 We show that by application of an applied magnetic field it might be possible to reorient the ferromagnetic 
moment through successive 90° rotations, which would then induce similar rotations of the spontaneous polarization. 
Likewise, application of an external electric field could reorient the spontaneous polarization which in turn would 
reorient the wave vector and thereby reorient the weak ferromagnetic moment. 

•5 Here we analyzed behavior near the phase transitions using mean field theory. However, non mean-field critical 
exponents can be accessed experimentally, as has been done for Ni^V^Os. [45j| 

• 6 We have given a detailed enumeration (see Fig. IT2|) of the domains arising from different realization of the 
OP's. We distinguish between domains whose bulk structure is macroscopically identifiable and those (similar to 
antiferromagnetic domains) that arise from a difference in phase that is not macroscopically accessible. We propose 
the association of domain walls with planes of "minimum contact" between octahedra. The question of the formation 
and dynamics of domains (especially those of phase domains) is broached. 
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Appendix A: Crystal Structure for Q^(qi) and Q^(qi) 

Here we verify that the crystal structure when the OP's for ireps and X^ at wave vector qi are simultaneously 
nonzero is Cmc2i. (This result also applies when the OP's are both at wave vector q2.) From Table [TT1 we have the 
characters listed in Table IIVI 

Now which operators transform like unity under both irreps? We see that we may choose 

TiT 2 , TiT^ 1 , Tf 1 T 3 , m z , m d Ti . 
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These indicate that we new primitive lattice vectors are 

ai = (1,1,0), a 2 = (l,T,0), a 3 = (-1/2, 1/2, 1/2). 

Also 

m d = (x,y,z,) -> (x,y,z) , 
m d Ti = (x,y,z) ->• (y,x+ l,z) 

TABLE IV: Characters of generators of the little group for irreps X% and X$ at wave vector qi. 



o = 


1 m d m z Ti T 2 T-s 




1-11-1 -1 -1 
-1 -1 1 -1 -1 -1 



To make contact with Ref. Ila we transform to orthorhombic coordinates: 



, x-y 1 , x+y 



In this coordinate system 



and the mirror operations are 



(0,0,1), ^ = (0,1,0), 4 = (1/2,-1/2,0) 



TO, 



(x',y',z') ->■ {x', y',z') 



[mrfTr]' = (x',y',z')^(x',y',z' + l/2) 

which coincides with the specification of space group Cmc2i in Ref. O One might object that we have not taken 
into account the fact that Eq. ([7]) indicates the presence of irrep . What that means is that this irrep is always 
allowed in Cmc2i. 

Appendix B: CRITICAL BEHAVIOR OF THE DIELECTRIC CONSTANT 

The following discussion parallels that given[38[ for the dielectric anomaly in Ni3V20s- We take the free energy (in 
the qi channel) for T near the lower transition [where Qa(qi) is already nonzero] to be 

T = y(T-T 2 )g 2 ( qi ) 2 + y(T-T 5 )[Q 5 , 1+ Q 5 , 2 ] 2 

WQ 2 (qi)[Q 5 ,i + Q S ,2] - V2X[P X Q 5<1 + P y Q 5 . 2 } 
+ lx E 1 P 2 - A / <Q 3 (qi))Q 2 (qi)[ J Px + P y ] 

-[P X E X + PyEy] . (Bl) 

The coupling terms proportional to A and A' suggest that the dielectric susceptibility will be singular at the lower 
transition where Q2 and Q5 appear. In this appendix we show this explicitly. We write 

PxQ$,l + PyQb,2 = 7: [Px + Py][Q 5,1 + Q 5,2] 



" Py}[Q5,i — Q5,2\ 



and 



PxPx PyPy 



[Px + Py][E x + E y 



1 



P x ~ Py][E x - Ey] 
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Note that in the qi channel Q$ t i — Q5.2 = and can be dropped. Also, since E x — E y and P x — P y do not couple to 
a critical variable we drop them too. Now set 



Thus the above free energy is 



P = [P x + P y ]/V2 , E=[E X + E v ]/V2 , 

Q5 = £?5,1 + Q5,2 ■ 



T = ~(T -T 2 )Ql + ^(T -T 5 )Q 2 5 +w'Q 2 Q 5 

where A" = A'(Q3(qi))-\/2- Now minimize with respect to Q2 and Q5 to get 

dF 

— = a'(T-T 2 )Q2+w'Q 5 -X"P = , 

OL} 2 

dT 

— = w'Q 2 + a"{T - T 5 )Q 5 - XP = , 
0V5 



so that 



where 



Q 2 = [\"a"(T-T 5 )-w'\]P/D 
Q 5 = [-A'V + a'(T - T 2 )\)P/D , 



D = a'(T-T 2 )a"(T-T 5 )-w' 2 
Then the equation for P is dF/dP = 0, or 

-\"q 2 -\q 5 + xe 1 p = e , 

so that the dielectric susceptibility is 

P 



XE = E 
= D 



'a"x E \T-T2){T-T 5 )-w' 2 X - E 1 



-X 2 a'(T - T 2 ) - X" 2 a"{T - T 5 ) + 2w' XX" 

Note that xe has poles at T — T 2 (which is close to T 2 ) and at T5 (which is close to T5). When the effects of w' , X, 
and A" can be treated perturbatively with respect to T 2 and x^J 1 , we find that 

T2 = T2 + ^— + ^—^^T2 + 6T2 
a a" ±2 a 

a result which is reasonable considering the couplings proportional to w' and A' in Eq. (Bl). For T near T2 we can 
therefore write 

~ = Dxe = Axe 

XE a'a"(T -f 2 )(T -f 5 ) T -f 2 ' 

where, in terms of D{T), we have 

^ D(T = f 2 ) _ a'a"T 2 (5T 2 ) - w' 2 _ X" 2 Xe 
a'a"T 2 a'a"T 2 a' ' 

which is a small amplitude attributable to the existence of the mixing term proportional to A' in Eq. (|B1|) . But this 
result does confirm the expected divergence in the dielectric constant at the lower transition where the polarization 
first appears. 
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